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ABSTRACT 

We present unambiguous evidence for a parsec scale wind in the Broad-Line 
Radio Galaxy (BLRG) 3C 382, the first radio-loud AGN whereby an outflow 
has been measured with X-ray grating spectroscopy. A 118 ks Chandra grating 
(HETG) observation of 3C 382 has revealed the presence of several high ionization 
absorption lines in the soft X-ray band, from Fe, Ne, Mg and Si. The absorption 
fines are blue-shifted with respect to the systemic velocity of 3C 382 by —840 ± 
60kms -1 and are resolved by Chandra with a velocity width of a — 340 ± 
70kms -1 . The outflow appears to originate from a single zone of gas of column 
density Nu — 1.3 x 10 21 cm -2 and ionization parameter log(£/ergcms -1 ) = 2.45. 

From the above measurements we calculate that the outflow is observed on parsec 
scales, within the likely range from 10 — 1000 pc, i.e., consistent with an origin in 
the Narrow Line Region. Finally we also discuss the possibility of a much faster 
(0.1c) outflow component, based on a blue-shifted iron Ka emission fine in the 
Suzaku observation of 3C 382, which could have an origin in an accretion disk 
wind. 

Subject headings: Galaxies: active — galaxies: individual (3C382) — X-rays: 
galaxies 
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1. Introduction 

At least 50% of radio-quiet AGN exhibit evidence for photoionized outflows in their 
X-ray spectra (Reynolds 1997; George et al. 1998; Crenshaw et al. 2003; Porquet et al. 
2004; Blustin et al. 2005; McKernan et al. 2007). The signatures of these winds consist 
of absorption and emission features at soft (0.4-2 keV) and hard (6-8 keV) X-ray energies, 
coincident with ionized O, N, Ne, Mg, Si and Fe lines blueshifted in the observer’s rest-frame. 
The inferred velocities of the winds are typically in the range 100 — 1000 km s -1 , but can 
be as high as ~ 0.1c in some sources (Chartas et al. 2002, 2003; Pounds et al. 2003; Reeves 
et al. 2003, 2008). It is also possible that the energy budget of the outflows of some AGN 
can approach a significant fraction of the bolometric or even Eddington luminosity (King &; 
Pounds 2003). 

In stark contrast, the X-ray evidence for nuclear outflows is very scarce in Broad-Line 
Radio Galaxies (BLRGs) and in radio-loud AGN generally. Previously, the quasar 4C +74.26 
showed weak absorption features at ~ 1 keV with ASCA (Ballantyne 2005), while two radio 
galaxies, 3C 445 and 3C 33 (Sambruna et al. 2007; Evans et al. 2006) also exhibit soft X-ray 
emission lines below 2 keV, which could originate from spatially extended material (in 3C 33; 
Torresi et al. (2009a)). Disk winds are also expected in radio-loud AGN as ingredients for 
jet formation (Blandford &; Payne 1982). 

In this Letter we present direct evidence for outflowing gas from the nucleus of the 
nearby (z = 0.05787), bright BLRG 3C 382. A re-analysis of our 118 ks Chandra HETG 
(High Energy Transmission Grating) observations (Gliozzi et al. 2007) revealed several blue- 
shifted absorption lines between 0.7 — 2.0 keV which suggest the presence of a large-scale 
(10 — 1000 pc) outflow in this source with a velocity of 800 km s -1 . Moreover, our 100 ks 
Suzaku spectrum of 3C 382 shows possible evidence for a blue-shifted Fe K emission line at 
7.5 keV (object rest- frame), implying nuclear gas outflowing with velocities up to 0.1c. 

The organization of this Letter is as follows. In § 2 we describe the Chandra and Suzaku 
data reduction and analysis; in § 3 the results of the spectral analysis; Discussion and 
Conclusions follow in § 4. Throughout this paper, a concordance cosmology with Hq = 71 
km s -1 Mpc -1 , Q a = 0.73, and D m =0.27 (Spergel et al. 2003) is adopted. Errors are quoted 
to 90% confidence for 1 parameter of interest (i.e. Ay 2 or AC — 2.71). 
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2. The Data 
2.1. The Chandra HETG 

Chandra observed 3C382 with the HETG for a net exposure of 118 ks between 27-30 
November 2005. The ±1 order spectra were summed for the MEG (Medium Energy Grating) 
and HEG (High Energy Grating) respectively, along with their response files. The summed 
first order count rates for the MEG and HEG are 0.867 counts s -1 and 0.379 counts s -1 respec- 
tively, while the MEG data were fitted between 0.5-7.0keV and the HEG from 1.0-9.0keV. 


2.2. The Suzaku XIS 

Suzaku observed 3C 382 between 27 and 30 April 2007, with a total net exposure with 
the X-ray Imaging Spectrometer (XIS, Koyama et al. (2007)) of 116 ks. The observation was 
taken in the XIS nominal pointing position. The revision 2 pipeline data were used with 
standard screening, while source spectra were accumulated from a 2.9' radius circle centered 
on the AGN, with a background spectrum taken from a region offset from the source. The 
front illuminated XIS 0 and XIS 3 were found to be consistent within 5% cross-normalization 
uncertainties and were subsequently combined, along with the response files, to maximize 
the signal to noise ratio. The total count rate of the co-added XIS 0+3 (hereafter XIS FI) 
spectrum is 5.10 counts s -1 . The spectrum was re-binned to 28.8 eV per energy channel. 
Data from the back illuminated XIS 1 were not considered here, as the spectra have lower 
signal in the Fe K band, but it was found to be consistent with the XIS FI. An analysis of 
the data taken from the Suzaku HXD (Hard X-ray Detector, Takahashi et al. (2007)) will 
be presented in a subsequent paper (Sambruna et al. 2009, in preparation). 


3. Results 

3.1. The warm absorber 

The Chandra HETG data were first analyzed by Gliozzi et al. (2007), who focused on 
the continuum and its variability. Our results are in agreement with theirs. Specifically, the 
MEG and HEG data were fitted by an absorbed power law with photon index T = 1.66+0.01 
plus a blackbody with kT — 92 ± 6 eV to parameterize the soft excess below 1 keV (Gliozzi 
et al. 2007), absorbed by a Galactic line of sight column of An,Gai = 7.0 x 10 2O cm -2 (Dickey 
&; Lockman 1990). Figure 1 shows the broad-band HETG spectrum fitted with an absorbed 
power-law only, to illustrate that a soft excess is clearly present below IkeV. The 0.5-9 keV 
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band flux is 6.4 x 10 -11 erg cm -2 s _1 . Even upon adding a blackbody to parameterize the soft 
excess, the fit is still formally unacceptable (y 2 /dof = 632/477, null probability 2.4 x 10 -6 ) 
as there are clear residuals around 1 keV that indicate the presence of a warm absorber. 

To analyse the warm absorber in detail, the HEG and MEG spectra were binned more 
finely to sample the resolution of the detector, at approximately HWHM the spectral resolu- 
tion (e.g. AA = 0.01 A bins for the MEG). For the fits, the C-statistic was used (Cash 1979), 
as there are fewer than 20 counts per resolution bin. The absorption lines were modelled 
with Gaussians and the continuum model was adopted from above. Table 1 lists the detected 
lines with their observed and inferred properties, and their significance as per the C-statistic. 
Figure 2 shows the portions of the HETG spectrum containing the strongest lines, with the 
model overlaid. 

The seven absorption lines in Table 1 and Figure 2 are all detected at high confidence 
(corresponding to AC > 18, or > 99.9% confidence for 2 parameters of interest). The lines 
likely arise from the Is — 2 p transitions of Neix, Nex, Sixm, and Sixiv and the 2 p — 3 d 
lines of Fexix-xxi. The two statistically weaker Is — 2 p lines of Mgxi and Mgxn may also 
be present, which have outflow velocities consistent with the other lines. 

Initially we assume that the lines have the same velocity width within the errors. The 
velocity width of the absorption lines is then a — 340±70kms _1 (or 780±160kms _1 FWHM) 
and the lines are clearly resolved. Even at 99% confidence (A C — 9.2 for 2 parameters), the 
velocity width is constrained to a — 340± 140 km s -1 . Upon allowing the velocity width of the 
individual lines to vary, then they are constrained to lie within the range a — 250—500 km s -1 
as shown in Table 1. The mean outflow velocity is — 810kms _1 . The overall fit statistic is 
C = 3804 for 3811 bins. 

We used the photoionization code XSTAR (Kallman et al. 2004) to derive the parameters 
of the absorber, assuming the baseline continuum described above, including the soft excess. 
Solar abundances are assumed throughout (Grevesse & Sauval 1998). An important input 
parameter is the turbulent velocity, which can effect the absorption line equivalent widths and 
hence the derived column density. We experimented with two different values of the turbulent 
velocity chosen to represent two likely extremes: (i) a lowest value of v tU rb — lOOkrns -1 
and (ii) u tU rb — 300kms _1 , the latter being consistent with the measured width of the 
absorption lines. The fitted continuum parameters are T = 1.68±0.02 and for the blackbody, 
kT — 110 ± 8eV. For case (i), then 7V H = (3.2 ± 0.6) x 10 21 cm -2 , the ionization parameter 
is log £*= 2.45^o;og and the outflow velocity is u ou t = — 8IO+55 kms -1 . The fit statistic is 
G/bins = 3795/3811. For case (ii), then N H = (1.30 ± 0.25) x 10 21 cm- 2 , log£ = 2.45lg;g? 


1 The units of £ are erg cm s 1 . 
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and the outflow velocity t> out = —840^0 kms -1 . The fit statistic is (7/bins = 3783/3811. If 
the warm absorber is not included in the model, then the fit statistic is substantially worse 
by AC = 220 (compared to model (ii)). Only a single outflowing layer of gas is required to 
model the warm absorber. 

The higher turbulence velocity model is statistically preferred and is consistent with 
the measured 340 kms -1 widths of the lines. Either model yields an outflow velocity of 
—800 kms -1 within a statistical error of < 10%. Hereafter we adopt the parameters from 
model (ii), as the turbulent velocity is consistent with widths of the individual absorption 
lines. 

However in neither model is the fitted column density of the absorber as high as the 
value of Nu ~ 3 x 10 22 cm -2 reported by Torresi et al. (2009b) from an analysis of a short 
34.5 ks XMM-Newton/RGS observation of 3C382 on April 28, 2008. If the column density 
of the warm absorber is fixed to the value of Ah = 3 x 10 22 cm -2 in the HETG spectrum, 
then the fit statistic is considerably worse (C /bins = 8390/3811). 

As a consistency check, we analyzed the archival RGS data of 3C 382 with model (ii) 
and the same continuum form as above. We found that the column density and ionization 
parameter were degenerate with each other, given the short exposure of the RGS spectrum. 
Thus for the RGS, A H = I.4+J3 x 10 22 cm -2 , log£ = 3.4 ± 1.0 and the outflow velocity 
u ou t — —1200^5oo kms -1 . The fit statistic improves only by AC — 25 to C/bins = 826/820 
upon adding the absorber. Thus within the larger errors, the parameters are consistent 
with those obtained from the HETG. The individual lines detected in the Chandra HETG 
observation were also compared to the absorption lines claimed by Torresi et al. (2009b) on 
the basis of the RGS data. With the higher statistics of the HETG compared to the RGS, 
we only confirm the detection of one line reported by Torresi et al. (2009b), the Fexx line 
at a rest-frame energy of 1025 eV 2 


3.2. The Blue-shifted Iron K line 

The Suzaku XIS FI data of 3C 382 were examined to measure any emission or absorption 
in the iron K band. The continuum model from the HETG analysis was adopted, while the 
warm absorber described above was also included in the analysis. A Compton reflection 
continuum produced via scattering off neutral matter also was included, using the pexrav 


2 The line reported at 1.356 keV by Torresi et al. (2009b) may also be associated with the Is — 2 p line of 
Mg xi, as noted in Table 1. 
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model (Magdziarz &; Zdziarski 1995), with a fixed reflection fraction of R — 0.5 3 for an 
inclination of cos 6 — 0.87 and an exponential cut-off of 200 keV. The reflection parameters 
were determined from an analysis of the Suzaku XIS, HXD and Swift BAT broad-band 
spectrum from 0.6-100 keV, the details of which will be presented in a subsequent paper 
(Sambruna et al. 2009, in preparation). The continuum flux measured by the Suzaku 
XIS was at a similar level to the Chandra observation, of 6.5 x 10 -11 erg cm -2 s -1 over the 
0.6-10 keV band. The power-law photon index was a little steeper than in the Chandra 
observation, with T = 1.85 ± 0.01, while for the blackbody component, kT — 84 ± 9eV. 

While the above model provides an adequate description of the continuum in the XIS 
band from 0.6-10 keV, clear residuals are observed near 6.0 keV and 7.1 keV, as shown in 
Figure 3, while the fit statistic is unacceptable (y 2 /dof = 490/306, null probability of 1.2 x 
10 - 1 °) These residuals correspond to 6.4 keV and 7.5 keV in the rest frame of 3C 382. The 
6.4 keV line is at a rest-frame energy of 6.41 ± 0.01 keV and has an equivalent width of 
48 ± 10 eV, while the line is resolved with a width of a = 98 ± 30 eV (corresponding to a 
FWHM of 10500 kms -1 ). The weaker 7.5keV line has a rest frame energy 7.51 ± 0.02keV, 
blueshifted with respect to the iron Ka band, with an equivalent width of 20 ± 7 eV and is 
unresolved (a < 0.1 keV). 

The fit statistic was significantly improved upon the addition of the two emission lines 
(y 2 /dof = 355/301, null probability of 1.7 x 10 -2 ). However the improvement in y 2 was 
smaller (Ay 2 = 21) upon adding the weaker 7.5 keV emission line. We therefore tested its 
statistical significance by performing 1000 Monte Carlo simulations under the null hypothesis 
of there being no blue-shifted line, with initial input continuum parameters as above. The 
method used is identical to that described in Markowitz et al. (2006) or Reeves et al. (2009) 
and we conduct a blind trial allowing for line emission at any arbitrary energy over the 
energy range 4-9 keV. Only 1 out of 1000 trial spectra show a Ay 2 as high as Ay 2 = 21 
hence the 7.5 keV emission is detected with a null probability of 1 x 10 -3 . 

The 6.4 keV line arises from near-neutral iron Ko; emission (i.e. Fe<xvil), with a 
velocity width of ~ 10000 km s -1 , implying an origin much closer than the warm absorbing 
gas (see below) of order 1000 gravitational radii. This is consistent with the outer accretion 
disk or Broad Line Region (BLR), noting that the Ho; FWHM for 3C382 is 11800 km s -1 
(Eracleous &; Halpern 1994). 

The origin of the 7.5 keV line, if real statistically speaking, is less clear. The closest likely 
atomic transition would be from neutral Ni Ko; at 7.45 keV, however given the expected Solar 


3 The reflection fraction R is defined as R = £l/2ir, where fl is the solid angle subtended by the Compton- 
thick matter to the X-ray source. Thus R = 1 for a plane parallel slab. 
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abundance ratio of Fe/Ni of 25-30, it is unlikely that Ni Ka would contribute significantly. A 
more likely identification would be with Fexxv or Fexxvi Is — 2 p at 6.70 keV and 6.97 keV 
respectively, although either case would require a blueshift of 0.08 — 0.12c. 


4. Discussion and Conclusions 

The Chandra HETG spectra have revealed an ionized outflow in the BLRG 3C382. 
The outflow parameters are well determined, with Ah — (1.30 ± 0.25) x 10 21 cm -2 , log £ = 
2.451 q;o 7 an d v out = —840^0 kms -1 , while the absorption line widths are resolved with 
a — 340 ± 70 km s -1 . 

To characterize the outflow, we define the (unabsorbed) ionizing luminosity L; on , which 
in XSTAR is defined from 1 — 1000 Rydberg. This depends on the continuum model fitted 
to the data and it is important to take into account the soft excess. Using the best fit 
powerlaw plus blackbody continuum, then L ion = 1.2 x 10 45 ergs -1 . We nore that if a 
different continuum form is used to parameterize the soft excess, e.g. a broken power-law, 
then the ionizing luminosity can be a factor of ~ 2 higher; keeping this caveat in mind we 
adopt the more conservative lower luminosity value of 1.2 x 10 45 ergs -1 


4.1. The Location of the Absorber 

The upper-bound to the wind radius (i? 0 ut) is determined by geometrical constraints, 
i.e. if the thickness of the absorber A R/R << 1 (valid for a thin shell) and as Ah = nAR, 
where n is the electron number density. As the ionization parameter of the absorber is defined 
as £ = L ion /nR 2 , then by substitution R out << T ion /A H £ = 3.3 x 10 21 cm (or << lkpc) 4 . 
The lower bound is formed by the escape velocity, i.e. for the gas to escape the system as 
an outflow then R esc > c 2 /v 2 R s , where R s = 2 GM/c 2 is the black hole Schwarzschild radius 
and v — 800 km s -1 . If for 3C382, M — lx 10 9 M 0 (with a 40% uncertainty, see Marchesini 
et al. (2004)), then R esc > 1.4 x 10 5 R S > 4.2 x 10 19 cm > 13pc. In other words the location 
of the soft X-ray outflow is likely bounded between 10 pc and 1 kpc. 

Furthermore the outflow velocity and absorption line FWHM of ~ 800 km s _1 is similar 
to the width of the narrow optical forbidden lines of [Om], [Ol] and [Sill], which for 3C382 
lie in the range 400 — 600kms _1 , as measured from the spectra of Eracleous & Halpern 
(1994). Thus the origin of the X-ray absorption appears consistent with the optical Narrow 


4 Note the same radius is obtained by integrating down the line of a sight of a homogeneous radial outflow. 
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Line Region (NLR). The coincidence between the soft X-ray absorbing gas and extended 
[O ill] emission has also been noted in the radio-quiet quasar, MR 2251-178 (Kaspi et al. 
2004). 

The density of the outflow is then n — L ion /£R 2 . For the parameters above, then 
n — 0.4 — 2400 cm -3 . While loosely constrained, this is consistent with typical expected 
NLR densities of ~ 10 3 cm -3 in AGN (Koski 1978). Note the absorber thickness is then in 
the range A R/R — 0.01 — 1. The mass outflow rate for a uniform spherical flow is M out — 
4irnR 2 m p v ou t , where nR 2 — Li on /£ and m p is the proton mass. Hence for the measured warm 
absorber parameters for 3C382, M out = 7.2 x 10 27 gs -1 or M out = 100 M 0 yr -1 . However 
this value may be considerably smaller if the gas covers < 4ir steradians solid angle or if it 
is clumped. 

Interestingly the mass outflow rate may be close to that measured in the high luminosity 
radio-quiet quasar PDS 456 (Reeves et al. 2009) which has a similar 10 9 M 0 black hole mass, 
except in PDS 456 the outflow velocity is closer to v = 0.2c. One possibility is that high 
velocity gas, launched from the innermost accretion disk is decelerated once it reaches pc to 
kpc scales. Intriguingly the putative blue-shifted Fe K emission detected in Suzaku could be 
evidence of a much faster 0.1c outflow in 3C382, observed on the scale of the accretion disk 
(i.e. R e sc ~ 100i? s .) Alternatively if very highly ionized outflowing gas exists, in order to 
produce blue-shifted Fe xxv or Fe xxvi emission in the Suzaku iron K band spectrum, then 
this gas may play a crucial role in collimating the jet (Konigl & Kartje 1994). 

The kinetic power of the outflow (for v oui = — 800 km s -1 ) is then E = 1/2 M out v 2 ut < 
2 x 10 43 ergs -1 , which energetically is a fairly insignificant 2% of the ionizing luminosity and 
is unlikely to contribute significantly towards AGN feedback (King 2003). However this may 
not be the case if there is indeed a 0.1c (disk-scale) wind closer in, as perhaps is suggested 
by the blueshifted iron K emission in the Suzaku data. 


This research has made use of data obtained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC), provided by NASA’s Goddard Space Flight Center. 
R.M.S. acknowledges support from NASA through the Suzaku and Chandra programs. We 
thank Tahir Yaqoob for assistance with the Chandra data analysis. 
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Fig. 1. — Chandra HETG MEG (black) and HEG (red) spectra of 3C382, binned coarsely 
at four times the resolution of the gratings, in order to show the broad-band continuum 
spectrum. The data is plotted as a ratio against an absorbed power-law, of photon index 
T — 1.66. A clear excess of counts is present in the soft X-ray band below 1 keV, while 
significant residuals below unity are also present between 0.7-1.0keV, indicating that a 
warm absorber is present. 
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Fig. 2. — Fluxed HETG spectra showing the comparison between the data binned at HWHM 
resolution and the best- fit absorption line model (solid line) described in the text. Several 
absorption lines are present in the HETG spectrum, as labelled in the above figure and listed 
in Table 1. 
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Fig. 3. — The Suzaku XIS FI spectrum plotted as a ratio to the best-fit continuum model 
described in the text, shown over the iron K band. A clear excess of emission is present at 
6.05 keV and 7.10 keV (observed frame), indicating emission lines at 6.40 keV and 7.50 keV 
in the rest-frame of 3C 382. The 6.4 keV line is the neutral emission from iron Ka, while the 
weaker 7.5 keV line might originate from Fexxv or Fexxvi, blueshifted by 0.08 — 0.12c. 


Table 1. Summary of HETG absorption line parameters. 
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a Observed energy of absorption line in eV. 

b Energy of absorption line in rest-frame of 3C382, in units eV. 

c ine identification and lab frame energy in eV in parenthesis. Atomic data are from 
http://physics.nist.gov 

Equivalent width, units eV. 

e Outflow velocity of absorption line, in units kms -1 . 
f lcr velocity width of absorption line, in units kms -1 . 
improvement in C-statistic, upon adding line to model. 
h Parameter is fixed in the model. 



